Number statistics of molecules formed from ultra-cold atoms 
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We calculate the number statistics of a single-mode molecular field excited by photoassociation or 
via a Feshbach resonance from an atomic Bose-Einstein condensate (BEC), a normal atomic Fermi 
gas and a Fermi system with pair correlations (BCS state). We find that the molecule formation 
from a BEC is a collective process that leads for short times to a coherent molecular state in 
the quantum optical sense. Atoms in a normal Fermi gas, on the other hand, are converted into 
molecules independently of each other and result for short times in a molecular state analogous to 
that of a classical chaotic light source. The BCS situation is intermediate between the two and goes 
from producing an incoherent to a coherent molecular field with increasing gap parameter. 

PACS numbers: 33.90.+h,03.75.-b,05.45.-a 
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A remarkable development following the availability of 
ultracold, quantum-degenerate atomic samples has been 
the coherent generation of molecular dimers by means 
of Feshbach resonances 0, 0> El or vlSl two-photon Ra- 
man transitions Qj from atomic Bose-Einstein conden- 
sates (BEC). These trailblazing experiments were soon 
followed by the generation of molecular dimers from two 
fermionic atoms 0, S 1 m which case the capability to 
tune the strength and change the sign of the two-body 
interaction between fermions opens up the possibility to 
carry out detailed studies of many longstanding ques- 
tions concerned with the so-called BEC-BCS crossover 
H, HI 0, 0| . The observation of Feshbach resonances 
between atoms of different species 0, also opens up 
the possibility of generating quantum-degenerate gases 
of heteronuclear molecules. These experiments have trig- 
gered a number of theoretical studies aimed at a better 
understanding of the process of coherent molecule for- 
mation and the properties of the coupled atom-molecule 
system 0, 0, ^ , but the description of the dynamics 
of these systems remains incomplete. In particular, the 
dependence of the quantum statistical properties of the 
resulting molecular field on those of the initial atoms is 
largely unknown. 

In this letter we apply methods of quantum optics to 
the coupled atom-molecule system and analyze the num- 
ber statistics of the molecules formed when starting from 
an atomic BEC, a normal Fermi gas, and a Fermi system 
with pairing (BCS type state). We restrict ourselves to 
zero temperature T — throughout, and further assume 
that the molecules are formed in a single mode, typically 
the ground state of a trap. The statistical properties of 
the molecular field are then determined by numerically 
integrating the relevant Schrodinger equations, and fur- 
ther analytical insight is provided by time-dependent per- 
turbation theory considerations. We find that the statis- 
tics of the molecular field provide a distinct signature of 
the initial atomic state, thereby offering a powerful diag- 
nostic tool to measure atomic statistical properties such 
as e.g. the presence of BCS pairing. 

Consider first a cloud of weakly interacting bosons well 
below the condensation temperature T c . It is a good ap- 
proximation to assume that all atoms are in the conden- 



sate, described by the condensate wave function Xo{%) 
with energy eigenvalue /i. 

Atom pairs can be transformed into trapped molecules 
by means of photo-association or by means of a Feshbach- 
resonance. The coupled system of atoms and molecules is 
described by the effective two-mode Hamiltonian fvA ITsj 
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where a) and c, & are the bosonic annihilation and cre- 
ation operators for the atoms in the condensate and for 
the molecules, respectively, 8 is the detuning between the 
molecular and atomic level, g is the effective coupling 
constant, and we set h = 1 throughout. 

Typical experiments start out with all the particles in 
the condensate and no molecules, corresponding to the 
initial state 



m = o)) = 
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where A^ a 

maximum possible number of molecules and |0) is the 
vacuum of both molecules and atoms. For simplicity we 
only consider the case of even numbers of atoms, but 
our treatment can easily be extended to odd atom num- 
bers. The Hamiltonian clearly conserves the total 
number of free and bound atoms, 2a) a + &c. Therefore, 
the evolution of the system from the initial state J2J) can 
be described in the basis 
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with Af being a normalization constant. The state \<j) n ) 
corresponds to n molecules and 2(iV max — n) atoms. 

Expanding the state of the system on this basis as 
I^W) = ^InVn^n) we can write the Schrodinger equa- 
tion as 
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dt 
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FIG. 1: Number statistics of molecules formed from a BEC 
with -/V max = 30 and 8 — 0. 



Solving this coupled set of N max equations numerically 
gives then the probability P n (t) = \y n (t)\ 2 of find- 
ing n molecules in the single-mode molecular field, the 
"molecule statistics." 

Figure shows P n (t) for 30 initial atom pairs and 
5 = 0. Starting in the state with zero molecules, a 
wave-packet-like structure forms and propagates in the 
direction of increasing n. Near N max the molecules be- 
gin to dissociate back into atom pairs. (Note that our 
two-mode model neglects "rogue photodissociation," an 
approximation appropriate for short enough times so that 
the condensate is not significantly depleted 01 •) 

We can gain further insight into the short-time dynam- 
ics of molecule formation by using first-order perturba- 
tion theory, which gives the mean molecule number 



n{t) = (fits) = (gt) 2 2N max (2N mlix - 1) 
and the second factorial moment 
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For Nxoax large enough we have g^ 2 \t\, t 2 ) — > 1, the value 
characteristic of a Glauber coherent field. From g^ and 
n(t) we also find the relative width of the molecule num- 
ber distribution as 



= Vg {2) + ' 



- 1. 



(7) 



It approaches n ' in the limit of large N max , typical 
of a Poisson distribution. This confirms that for short 
enough times, the molecular field is coherent in the sense 
of quantum optics. 

Let us now turn to the case of photoassociation 
from two different species of non-interacting ultra-cold 



fermions. The two species are denoted by spin up and 
down. At T = 0, the atoms fill a Fermi sea up to an en- 
ergy /i. Weak repulsive interactions give rise only to mi- 
nor quantitative modifications that we ignore in this pa- 
per. We refer to this system of non-interacting Fermions 
as a normal Fermi gas (NFG) [T^| . 

As before we assume that atom pairs are coupled only 
to a single mode of the molecular field which we assume to 
have zero momentum for simplicity. Only pairs of atoms 
with opposite momenta +k and —k and opposite spin can 
then be combined into molecules. Introducing pseudo- 
spin operators in the subspace of each pair of atoms with 
opposite momenta and spin as 
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1, 
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where Ck,H and c k ^ are the annihilation and creation 
operators for a fermion of momentum k with up and 
down spin, the Hamiltonian describing the coupled atom- 
molecule system takes the form |2(j, |21| 



NFG 



k \ k 



H.c. , (9) 



where Ek = k 2 /2M is the kinetic energy of an atom of 
momentum k and mass M and S is the detuning between 
the molecule energy and the Fermi energy. This Hamil- 
tonian is formally equivalent to the well studied Tavis- 
Cummings model of quantum optics, which describes the 
coupling of an ensemble of two-level atoms to a single 
mode of the electromagnetic field |22|. 

The atomic kinetic energies Ek, which in this anal- 
ogy correspond to the energies of the fictitious two-level 
atoms, give rise to inhomogeneous broadening. In case 
this broadening can be neglected the atoms are conve- 
niently described in terms of the eigenstates \S, m) of S 2 
and S z , 

S 2 \S,m) = S{S +l)\S,m), S z \S,m) = m\S,m), (10) 



where S 



l = x,y, z, and we have as usual 



*fk) = (%) + & (k)^ 2 and °{k) = ( & (k) - a (k)\/ 2i - In 
that limit the Hamiltonian © commutes with S z + a* a 
and, like in the BEC case, the system can be uniquely 
described by the probability amplitudes y n for finding n 
molecules. This degenerate model is studied in detail in 
Ref. [23|. We focus here instead on the situation where 
the kinetic energy of the atoms cannot be neglected. 

Figure [21 shows the molecule statistics obtained by a 
numerical integration of the Schrodinger equation corre- 
sponding to the Hamiltonian (JSJ. The result is clearly 
both qualitatively and quantitatively very different from 
the case of molecule formation from an atomic BEC. 

From the Tavis-Cummings model analogy we expect 
that for short times the statistics of the molecular field 
should be chaotic, or "thermal" , much like those of a 
single-mode chaotic light field. This is because each indi- 
vidual atom pair "emits" a molecule independently and 







FIG. 2: Number statistics of molecules formed from a normal 
Fermi gas. This simulation is for N a = 20 atoms, the detuning 
is 5 = and the Fermi energy is /i = O.lg. The inset shows fits 
of the number statistics to thermal distributions for various 
times. 



FIG. 3: Number statistics of molecules formed from a Fermi 
gas with pairing correlations. For this simulation the detuning 
is 8 = 0, the Fermi energy is /i = O.lg and the background 
scattering strength is V = 0.03g resulting in N a ~ 9.4 atoms 
and a gap of A ft; 0.15g. 



without any phase relation with other pairs. That this 
is the case is illustrated in the inset of Fig. [21 which 
fits the molecule statistics at selected short times with 
chaotic distributions of the form 



P, 



^/(n) 



n, thermal 



-n/(n) 
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The increasing 'pseudo-temperature' (n) corresponds to 
the growing average number of molecules as a function 
of time. 

We can again determine the short-time properties of 
the molecular field in first-order perturbation theory. For 
the mean number of molecules we find 



n(t) = (gt) 2 2N a . 



(12) 



It is proportional to N a , in contrast to the BEC re- 
sult, where n was proportional to N 2 , see Eq. Q. This 
is another manifestation of the independence of all the 
atom pairs from each other: While in the BEC case the 
molecule production is a collective effect with contribu- 
tions from all possible atom pairs adding constructively, 
there is no such collective enhancement in the case of 
Fermions. Each atom can pair up with only one other 
atom to form a molecule. 

For the second factorial moment we find 



gW(t u t 2 ) = 2[l 



1 



2N a 



(13) 



which is close to two, typical of a chaotic or thermal field. 

Unlike repulsive interactions, attractive interactions 
between fermions have a profound impact on molecule 
formation. It is known that such interactions give rise to 
a Cooper instability that leads to pairing and drastically 
changes the qualitative properties of the atomic system. 



The BCS reduced Hamiltonian including interactions be- 
tween the atoms is 241 



Hbcs = -Hnfg — &k ^k" 



(14) 
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The BCS ground state is found by minimizing (Hbcs 
liN) using the ansatz 



|BCS) = nK + ^ fc + )D, 



with the result 



(15) 

(16) 
(17) 



where = (Ek — (J.) is the kinetic energy of the atoms 
measured from the Fermi surface and A = V^2 k ui-vi- 
is the gap parameter. It is determined by numerically 
solving the gap equation 




V ^ 



1 



(18) 



which is readily done for the small atom numbers at hand. 
The dynamics is then obtained by a numerical integra- 
tion of the Schrodinger equation with |BCS) as the initial 
atomic state and the molecular field in the vacuum state. 

Figure [31 shows the resulting molecule statistics for 
V = 0.03g, which corresponds to A = 0.l5g for the sys- 
tem at hand. (This large background scattering strength 
was chosen in order for the gap equation to have a posi- 
tive solution for the small particle numbers to which we 
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FIG. 4: gr J (0 , ) as a function of the gap parameter A. 



are limited by computer memory requirements.) Clearly, 
the molecule production is much more effective than in 
the case of a normal Fermi gas. The molecules are pro- 
duced at a higher rate and the maximum number of 
molecules is larger. The evolution of the number statis- 
tics is reminiscent of the BEC case. 

The short-time dynamics is again obtained in first- 
order perturbation theory, which gives now 



n(t) = (gt) 2 [ ^ u k v k u k ,v k > 

2 



E«2 



(gt) 2 



(19) 
(20) 



In addition to the term proportional to N a represent- 
ing the incoherent contribution from the individual atom 
pairs that was already present in the normal Fermi gas, 
there is now an additional contribution proportional to 
(A/V) 2 . Since (A/V) can be interpreted as the number 
of Cooper pairs in the quantum-degenerate Fermi gas, 
this term therefore can be understood as resulting from 
the coherent conversion of Cooper pairs into molecules 



in a collective fashion similar to the BEC case. The co- 
herent contribution results naturally from the nonlinear 
coupling of the atomic field to the molecular field. This 
nonlinear coupling links higher-order correlations of the 
molecular field to lower-order correlations of the atomic 
field. For the parameters of Fig. [3] A/V « 6.5 so that 
the coherent contribution from the Cooper pairs clearly 
dominates over the incoherent contribution from the un- 
paired fermions. Note that no signature of that term 
can be found in the momentum distribution of the atoms 
themselves. Their momentum distribution is given by 
i^k a^k,a) = v \i a = T;I an d is very similar to that of a 
normal Fermi gas. The short-time value of g^(t\,t2), 
shown in Fig. 0] decreases from the value (|13[) for a nor- 
mal Fermi gas at A = down to one as A increases, 
underlining the transition from incoherent to coherent 
molecule production. 

In summary, we have applied concepts of quantum op- 
tical coherence to characterize the coherent generation 
of a molecular field created from a quantum-degenerate 
atomic sample. For atoms initially in a BEC the result- 
ing molecular field is to a good approximation coherent. 
This is in sharp contrast to the case of atoms in a normal 
Fermi gas, where we have made use of an analogy with 
the Tavis-Cummings model to show that the statistics 
of the resulting molecular field is very similar to those 
of a single-mode chaotic light field. The BCS case inter- 
polates between the two extremes, with an 'incoherent' 
contribution from unpaired atoms superposed to a 'co- 
herent' contribution from atomic Cooper pairs. We see, 
then, that the quantum statistics of the atomic sample 
has a profound impact on the quantum statistical prop- 
erties of the resulting molecular field. Conversely, these 
statistics provide a distinct signature of the initial atomic 
state and suggest the use of single molecule counting as 
a diagnostic tool for atomic states. 
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